Morphological Thermodynamics of Fluids: Shape Dependence of Free Energies 
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We examine the dependence of a thermodynamic potential of a fluid on the geometry of its 
container. If motion invariance, continuity, and additivity of the potential are fulfilled, only four 
morphometric measures are needed to describe fully the influence of an arbitrarily shaped con- 
tainer on the fluid. These three constraints can be understood as a more precise definition for the 
conventional term "extensive" and have as a consequence that the surface tension and other ther- 
modynamic quantities contain, beside a constant term, only contributions linear in the mean and 
Gaussian curvature of the container and not an infinite number of curvatures as generally assumed 
before. We verify this numerically in the entropic system of hard spheres bounded by a curved wall. 
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Although thermodynamics is built on extremely gen- 
eral assumptions its implications are far reaching and 
powerful. One basic building block is geometry which has 
a long history in thermodynamics and statistical physics 
of condensed matter. The formulation of thermodynam- 
ics in terms of differential forms [1] , scaled-particle theory 
(SPT) for fluids [2], depletion forces [3] of colloids in bio- 
logical cells [4], and density functional theory [5] (DFT) 
based on fundamental geometric measures [6-8] are only 
a few examples for the importance of a general geometric 
point of view on thermodynamic properties. In Refs. [9- 
11] the structure and phase behavior of microemulsions 
was explained assuming that the free energy of the fluid 
on a mesoscopic scale is given solely by four fundamen- 
tal geometric measures [12]. Here, we show numerically 
for the first time that these four fundamental measures 
are sufficient to describe accurately the free energy of 
a hard-sphere fluid in contact with a complexly shaped 
wall which supports the more general assumptions made 
in Refs. [9-12]. 

The grand potential Vl = £l[S;T, fi\ of a fluid depends 
on the temperature T and the chemical potential fi of 
the system, as well as on certain geometrical quantities 
which describe the shape of the container that bounds 
the system S. What are these thermodynamically rele- 
vant morphological parameters? One usually argues that 
every thermodynamic potential is an extensive quantity, 
which means that it scales linearly with the " size" of the 
system S. By partitioning a large system into identical 
smaller subsystems one normally assumes that Q[S;T, fj] 
is proportional to the volume V = V[S] of the system and 
uses as ansatz £l[S; T, fj] = w(T, /i) • V[S] . The intensive 
quantity w(T, /z) is the negative of the pressure p(T,fi) 
and is independent of the "size" of the confining con- 
tainer of S. This simple ansatz however is only valid for 
infinite bulk, i.e. "border-less" systems. If S is bounded 
by a container, f2 depends on the shape of the container 
in a potentially complicated manner and is convention- 
ally described by an infinite expansion in powers of the 
curvatures of the wall. However, we show that general 
considerations restrict this functional dependence on the 
shape to a linear combination of only four morphological 



measures, if all intrinsic (correlation) length scales are 
small compared to the system "size". This finding is in 
particular important for systems such as porous media 
[13], biological cells [4], or complex fluids such as mi- 
croemulsions [9] where fluids are confined by complexly 
shaped compartments and where the dependence of ther- 
modynamic quantities and transport properties on the 
shape of pores or cells has significant functional and bi- 
ological consequences. Our arguments can however not 
be applied, to critical phenomena, or if long ranged fluid- 
fluid or fluid- wall interactions are considered or if wetting 
or drying phenomena [14] occur at the wall, as intrinsic 
lengths in such systems have a macroscopic size. 

Now we focus on the dependence of Q[S] on the shape 
of the system S. The actual form of this mapping from 
a container onto a real number is given by the type and 
state of the fluid under consideration and is a compli- 
cated integral over the phase-space of the system, which 
usually can only be calculated approximately. However, 
we impose the following three physical restrictions on this 
mapping: 

(i) Motion invariance: The thermodynamic poten- 
tial of a system must be independent of its location and 
orientation in space, i.e. fi^S] = Sl[S] for all translations 
and rotations g in 3 dimensions. 

(ii) Continuity: If a sequence of convex sets S n 
converges towards the convex set S for n — > oo, then 
^[<Sn] - > £l[S]. Intuitively, this continuity property ex- 
presses the fact that an approximation of a convex do- 
main by e.g. convex polyhedra also yields an approxi- 
mation of the thermodynamic potential Q[S] by n[SVi]. 
In SPT it is shown that continuity is violated if the con- 
tainer S is very small and comparable in size to that of 
fluid particles [15]. Here we consider a larger container 
S. 

(iii) Additivity: The functional of the union Si U S2 
of two domains Si, i — 1,2, is the sum of the func- 
tional of the single domains subtracted by the intersec- 
tion: n[SiUS 2 ] = ^[5i]+fi[S , 2 ]-rj[S'inS'2]. This relation 
generalizes the common rule for the addition of an exten- 
sive quantity for two disjunct domains Si fl £2 = to the 
case of overlapping domains by subtracting the value of 
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the thermodynamic quantity of the double-counted inter- 
section. Note that the intersection Si n S2 does not need 
to be a volume but can rather be an area or a line for ad- 
jacent containers Sj. Additivity can break down if long 
ranged interactions are present or if the system develops a 
macroscopic intrinsic length scale. A fluid can be consid- 
ered additive even inside a concave container, if opposing 
walls are separated by several correlation lengths. 

Naturally the question arises about the most general 
form of a potential that satisfies these three conditions. 
The Hadwiger theorem [16, 17] states that every motion- 
invariant, continuous and additive functional in 3 dimen- 
sions can be written as a linear combination of the volume 
V = j s dV, the surface area A — J as dA, the integrated 
mean curvature C — \ ac HdA, and the Euler character- 
Therefore we write 



istic X = J gs KdA of the container. 

n[S] = -pV[S] + <tA[S] + kC[S] + RX[S] (1) 

as a complete expression for the grand canonical poten- 
tial, if the aforementioned conditions are satisfied [9- 
12]. The pressure p(T,p), the surface tension at the 
planar wall a(T,p), and the bending rigidities n{T, p) 
and k(T, p) are properties of the fluid and the wall-fluid 
interaction, but are independent of the actual shape of 
the bounding wall. The latter two thermodynamic co- 
efficients describe the influence of the curvature of the 
wall. Similar coefficients are also used for the Hclfrich 
Hamiltonian [18], which describes the free energy cost 
of bending a membrane. It is compatible with Eq. (1) 
on a length scale larger than the persistence length of 
the membrane, where renormalized contributions propor- 
tional to H 2 vanish due to thermal fluctuations. 

Note that Eq. (1) can be easily applied even to com- 
plexly shaped objects because the shape of S enters f2 
only via the four simple morphometric measures V , A, 
C, and X, while the thermodynamic coefficients a, k and 
R can be determined in a simple geometry. 

Thermodynamic quantities can be derived directly 
from SI and inherit a simple dependence on the shape 
of S by virtue of Eq. (1). The interfacial tension 7 = 
(Cl+pV) /A, which measures the total change in the grand 
potential per unit area introduced by the wall, can be 
evaluated using 



7 = a + kH + kK, 



(2) 



where H = C/A and K = X/A are the averaged mean 
and Gaussian curvatures of the bounding wall. These ge- 
ometrical quantities can be calculated from the principal 
radii of curvature R\ and R2 via H = (1/Ri + l/i?2)/2 
and K = l/(i?ii?2)- Note that this further justifies the 
ansatz used in SPT [2, 19] for the interfacial tension and 
shows that the analytic dependence of 7 on the curvature 
is a direct consequence of the additivity of the grand po- 
tential. No higher powers or derivatives of H or K con- 
tributing either to 7 or to ft. 

Closely related to the interfacial tension is T, the excess 
(over the bulk) amount of fluid adsorbed at the wall per 
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FIG. 1: The expansion coefficients as defined in Eq. (1) of 
the grand potential of a hard sphere fluid. For each value of 
r), the four thermodynamic coefficients can be used to calcu- 
late thermodynamic quantities for arbitrarily shaped systems. 
The solid line of the inset shows the relative error of 7 at a 
cylinder with radius R c calculated using the thermodynamic 
coefficients at ij = 0.3. This error can be compared to the 
numerical relative error for a sum rule [19, 20] (dashed line), 
which gives an estimate for the accuracy of our DFT data. 
Both errors are of the same order of magnitude such that our 
numerical data is in agreement with the prediction of Eq. (1). 
The fluid was modeled via Rosenfeld's FMT [6]. 



unit area, which is defined as T = (1/A) J s (p(r) — p) dV, 
where p is the bulk density and p(r) is the inhomogc- 
neous density distribution of the fluid. It inherits the 
morphometric form of the interfacial tension via Gibbs' 
adsorption theorem 
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If a fluid is bounded by a container one denotes the 
number density closest to the wall as contact density p c , 
which becomes p c , when averaged over the boundary sur- 
face. For hard walls p c can be regarded as a thermody- 
namic quantity because of an exact sum-rule [19, 20] and 
also features a simple dependence on the geometry: 



p c = p + 2aH + kK. 



(4) 



This relation can be derived by generalizing the argu- 
ments given by Henderson [20]. Note that the density 
distribution of the fluid away from contact depends in a 
more complicated way on the curvature of the wall. 

We test the vanishing of higher powers of H and K in 
thermodynamic quantities such as O, 7, T, or p c by exten- 
sive numerical studies of a fluid of hard spheres of radius 
R bounded by a hard wall using DFT. For a given chem- 
ical potential p, the bulk number density p of spheres or 
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FIG. 2: Curvature expansion coefficients for the contact den- 
sity of a hard-sphere fluid as a function of the packing frac- 
tion rj. Only the additive contributions p c P , p c H , and p c K con- 
tribute. The inset shows the relative error for the contact 
density at a spherical container with radius 57? between a di- 
rect DFT calculation and two approaches: For the solid curve 
all shown curvature expansion coefficients have been used, for 
the dashed one only additive contributions. The error, which 
is due to the numerical inaccuracies, is in both cases roughly 
the same. This shows that non-additive terms do not con- 
tribute to the thermodynamic quantity p c . Similar results 
were found for 7 (see Fig. 1). 



equivalently the packing fraction rj = 4irR 3 p/3 are fixed. 
Hard-core interactions do not introduce an energy scale 
such that we can account for the temperature dependence 
of all thermodynamic quantities by scaling energies with 
(3 = 1/(&bT). One can expect the grand potential of 
a hard-sphere fluid to be additive as hard core interac- 
tions are short ranged and the correlation length remains 
small in the fluid phase. Unfortunately, for arbitrarily 
shaped walls, there exists no reliable and accurate direct 
method to calculate thermodynamic quantities. How- 
ever, for simple geometries such as planar, spherical, or 
cylindrical walls, it is possible to apply e.g. DFT tech- 
niques [5], which allow the calculation of a, k and k. For 
these geometries, the curvatures and the contact den- 
sity are constants over the surface such that we can re- 
place averages over the surface by their local quantities. 
Rosenfeld's fundamental measure theory (FMT) [6] and 
the White-Bear version of FMT [7, 8] have proven to de- 
scribe the thermodynamics and structure of hard-sphere 
fluids very accurately and both give qualitatively equiv- 
alent results in the following analysis, as well as the non 
FMT-based Tarazona Mark I functional [21]. It is impor- 
tant to realize that minimizing a FMT functional does 
not restrict the results to standard SPT behavior [19]. 

The thermodynamic coefficients of Eq. (1) are shown 
in Fig. 1 as functions of the packing fraction 77 of the fluid. 



As dividing interface, which determines the surface area 
A as well as the curvatures H and K, we choose the sur- 
face where the contact density is measured. From these 
coefficients one can obtain values for 7 and p c for various 
geometries. In the inset of Fig. 1 we show the relative 
error (full line) of 7 of a hard-sphere fluid with packing 
fraction rj = 0.3 at a cylinder with radius R c as calcu- 
lated with the thermodynamic coefficients compared to 
that obtained directly from DFT. This relative error is 
of the same order of magnitude as the relative numerical 
error of the contact sum-rules, Eq. (4), indicating clearly 
that the very small deviation between the morphometric 
interfacial tension and that from a direct DFT calcula- 
tion is a numerical error of our calculation. 

In order to show that the linear expansion in H and K 
of 7 and p c is sufficient, we include for the analysis of our 
numerical data also non-additive terms and show that 
they do not contribute cither to 7 or to p c . To this end 
we introduce the curvature expansion of e.g. the contact 
density 

P c = Pp + PhH + P C K K + p^H 2 + p c HK HK (5) 

The coefficients p\ for £ = P,H, K, H 2 , HK, H 3 ,... can 
be determined numerically using DFT. The results in 
Fig. 2 show that irrespective of the bulk density only 
the additive terms p c Pl p c Hl and p c K contribute to p c . 
Note that in particular p c K (additive) contributes to the 
contact density whereas p c H2 (non-additive) does not al- 
though both K and H 2 feature the same quadratic de- 
pendence on the radii of curvature. This shows that the 
expansion in H and K is complete and not a truncated 
power series. In the inset of Fig. 2 we show the relative 
error of the contact density of a hard-sphere fluid at a 
sphere with radius 5R as a function of rj obtained by the 
morphometric form (dashed line) of p c , Eq. (4), and by 
the generalized curvature expansion (full line), Eq. (5), 
both compared to the direct calculation of DFT. Since 
in both cases the relative error is roughly the same we 
conclude that additive contributions are sufficient to de- 
scribe the influence of curvature on the contact density 
P c - 

To further highlight that the remarkably simple func- 
tional dependence on only four morphometric measures 
of the container is a peculiar feature of thermodynamic 
quantities, we introduce, in addition to the hard-wall po- 
tential, a short-ranged soft part to the wall-fluid poten- 
tial of the form V so f t {u) — Vbexp(— u/X), where u de- 
notes the normal distance from the dividing interface. 
It can be either attractive, if Vq < 0, or repulsive, if 
Vq > 0. Note that this additional part of the wall-fluid 
interaction leaves the morphometric form of thermody- 
namic quantities such as the grand potential f2, Eq. (1), 
and the surface tension 7, Eq. (2), unchanged. The con- 
tact density p c however loses its morphometric form for 
any non-zero value of Vo because the sum-rule Eq. (4) 
acquires an additional term [20] . 



4 




-1 -0.5 0.5 1 1.5 2 

PV 

FIG. 3: If the external potential of a hard container is 
perturbed by an additional potential as e.g. V so ft(u) — 
Vo exp(— u/X) the contact density p c is no longer a thermody- 
namic quantity for Vo 7^ 0. In this case non-additive contri- 
butions such as p c H 2 are necessary to describe the dependence 
on the curvatures of the wall. Here we show p c K and p c H2 in a 
logarithmic plot for A = R and a packing fraction r\ = 0.3. It 
is striking that even very small perturbations lead to p c H 2 / 
which is a clear indication that additivity of the contact den- 
sity p c is destroyed. In contrast, p c K contributes to the contact 
density for all Vo- 

This effect can be demonstrated by performing a cur- 
vature expansion, Eq. (5), of the contact density as a 
function of Vq. In the case of a purely hard wall, Vq = 0, 
the curvature expansion contains only additive contribu- 
tions (pp, p%, and p c K ), whereas for non-zero values of 
Vo this is no longer the case and the amplitude of e.g. 
p c H 2, increases with an increasing value of | Vol- As an 
example of this behavior we display in Fig. 3 p c K (addi- 
tive) and p c H 2 (non-additive) for a hard-sphere fluid with 
packing fraction -q = 0.3 for values of f3Vo from —1 to 2 
and A = R. It is very striking that already for very small 
amplitudes |Vb| 7^ 0, for which the soft part of the wall- 
fluid interaction can be regarded as a small perturbation, 
the contact density p c acquires additional non-additive 



contributions. We verified that the interfacial tension 7 
keeps its morphometric form (not shown), as expected. 

Monte-Carlo simulation confirm that the contact den- 
sity of a hard-sphere fluid with 77 = 0.3314 around a 
biaxial ellipsoid with half axes (4, 4, 10)i? as obtained by 
Eq. (4). Within the statistical errors of the simulation 
(about 1%), the results of Eq. (4) together with the ther- 
modynamic coefficients shown in Fig. 1 agree with the 
simulation data [22] . 

Furthermore, we also confirmed that Eq. (2) holds ac- 
curately in concave geometries if the perturbations intro- 
duced by the container do not interfere at caustic points. 
For a relatively high packing fraction of hard spheres of 
rj = 0.4 the deviation between Eq. (2) and direct DFT 
results is smaller than 1.5% for all radii of the cylinder 
larger than 5R. 

We presented a more precise definition than the con- 
ventional term "extensive" to describe the dependence 
of a thermodynamic quantity of a fluid with short in- 
trinsic length scales on the shape of the system. The 
assumption that the grand potential of a fluid is motion- 
invariant, continuous and additive allows an expansion 
in terms of only four simple morphological functionals 
[9-12]. As a consequence, the curvature expansion for 
thermodynamic quantities terminates after linear terms 
in mean and Gaussian curvature. This observation allows 
a calculation of thermodynamic quantities for complexly 
shaped objects with a greatly reduced effort in compar- 
ison to direct methods. The ideas presented here and 
tested numerically for the hard-sphere fluid should be 
applicable also to fluids with short-ranged interactions 
provided that intrinsic length scales are small compared 
to typical features of the container. If however, length 
scales are comparable, continuity or additivity may be 
violated and then thermodynamic quantities acquire ad- 
ditional contributions. 

The present approach is currently being extended to 
include local structure (density profiles) which will allow 
the prediction of entropic contributions to the effective 
interactions of complicated macromoleculcs. 

We thank H. Wagner, R. Evans and J.R. Henderson 
for stimulating discussions. 
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